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Abstract 

We theoretically study the resonant tunneling and Fano resonance in quantum dots with electron- 
phonon (e-ph) interaction. We examine the bias-voltage (V) dependence of the decoherence, using 
Keldysh Green function method and perturbation with respect to the e-ph interaction. With optical 
phonons of energy ujq, only the elastic process takes place when eV < Uq, in which electrons emit 
and absorb phonons virtually. The process suppresses the resonant amplitude. When eV > loq, 
the inelastic process is possible which is accompanied by real emission of phonons. It results in 
the dephasing and broadens the resonant width. The bias-voltage dependence of the decoherence 
cannot be obtained by the canonical transformation method to consider the e-ph interaction if its 
effect on the tunnel coupling is neglected. With acoustic phonons, the asymmetric shape of the 
Fano resonance grows like a symmetric one as the bias voltage increases, in qualitative accordance 
with experimental results. 

PACS numbers: 71.38.-k, 73.21.La, 73.23.-b 
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I. INTRODUCTION 



In semiconductor quantum dots, preservation of quantum coherence is an important issue 
for the application to the quantum information processing. To examine the coherence, the 
transport measurements have been reported using an Aharonov-Bohm (AB) ring with an 
embedded quantum dot, as an interferometer The AB oscillation of the current has been 
observed as a function of magnetic flux penetrating the ring. The conductance usually shows 
a Breit-Wigner resonance of symmetric Lorentzian shape, as a function of gate voltage which 
changes the electrostatic potential in the quantum dot.^ This indicates that electrons pass 
by a discrete level in the quantum dot coherently by the resonant tunneling. The transmitted 
wave of the electrons interferes with the wave traversing the other arm of the ring (reference 
arm), which results in the AB oscillation. 4 

When the higher-order interference between the waves is important, the resonant shape 
is modified. Kobayashi et al. have observed an asymmetric peak of the conductance, as a 
function of the gate voltage. 3 This is ascribable to the Fano resonance, which is generally 
caused by the interference between a discrete level and the continuum of states.— The Fano 
resonance is observable when the wave through a discrete level interacts with the wave 
traversing the reference arm sufficiently before it goes out of the ring: A naive picture in 
the path integral formalism is as follows. An incident wave from the source is split into two 
arms, one includes a quantum dot and the other is the reference arm, and meets at the other 
side of the ring. Some part of the wave goes out to the drain and the rest goes back to the 
original point along both arms and interacts, and so forth. If the phase relaxation length, 
1$, is comparable with the size of the ring, the shortest paths dominantly contribute to the 
current and result in the AB oscillation. Then we observe the Breit-Wigner resonance, as 
a function of gate voltage, which reflects the amplitude of a wave through the quantum 
dot. If is much larger than the ring size, the contribution from many paths modifies the 
resonant shape to the Fano resonance; an asymmetric shape with peak and dip is produced 
by positive or negative interference on either side of the Breit-Wigner resonance where the 
phase shift differs by tt from each other. In Ref. 0, the asymmetric shape of the current 
becomes a symmetric one when the bias voltage increases. This implies that the finite bias 
significantly reduces 1$. 

The resonances in the AB interferometer have been theoretically studied by several 
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groups . 6 i 7 i 8 i 9 i 10 i 1:L i 12 i 13 i 14 i 15 i 16 The electron-electron interaction has been investigated, which 
leads to the "Fano-Kondo effect."^ - Few theoretical work, however, has been done for the 
decoherence in the transport under a finite bias. In our previous paper, we have consid- 
ered the dephasing effect in the reference arm phenomenologically on the Fano resonance.— 
The dephasing effect on the AB oscillation by the spin-flip has been studied (in the case 
of Breit-Wigner resonance) when an odd number of electrons are localized in the quantum 
dot.^ 

In the present paper, we examine the decoherence in the resonances by the electron- 
phonon (e-ph) interaction. Although the e-ph interaction in quantum dots has been studied 
by several groups , 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 the bias- voltage dependence has not been clarified even 
in the Breit-Wigner resonance. Hence we discuss the e-ph interaction under a finite bias in 
the Breit-Wigner resonance as well as in the Fano resonance. The aims and main results of 
our study are as follows. 

(i) First of all, we study the effects of longitudinal optical phonons on the Breit-Wigner 
resonance. To calculate the current under a finite bias, we adopt the Keldysh Green function 
metho d 26127 and consider the e-ph interaction by the perturbation expansion . 17 ^ 8 Within the 
second-order perturbation, we distinguish elastic and inelastic processes at zero temperature. 
When the bias voltage is smaller than the phonon energy c^o, only the elastic process takes 
place, in which electrons virtually emit and absorb phonons. This process does not cause 
the dephasing, but it suppresses the resonant peak height without any change of the peak 
width. When the bias voltage is larger than uq, an inelastic process is possible, in which 
electrons actually emit phonons. This results in the dephasing 2 ^ and broadens the peak 
width. 

(ii) In some papers , 19 i 20 i 23 i 28 the e-ph interaction is treated by the canonical transformation 
method, in which e-ph interaction is considered exactly in the quantum dot, whereas its effect 
on the tunnel coupling is disregarded. We calculate the Breit-Wigner and Fano resonances 
with optical phonons by the method. By comparison with the perturbation calculation, we 
find some problems with the canonical transformation method to obtain the bias-voltage 
dependence. 

(iii) We examine the effects of acoustic phonons on the Breit-Wigner resonance and Fano 
resonance. For the purpose, we extend the self-consistent Born approximation to the finite- 
bias transport in the Keldysh Green function formalism . 17 > 18 We find strong suppression of 
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the peak height in the Breit-Wigner resonance since both elastic and inelastic processes are 
possible even at low bias in the case of acoustic phonons. In the Fano resonance, we show 
that the asymmetric shape is significantly influenced so that the dip almost disappears. This 
is in qualitative agreement with the experimental results.— 

Concerning (i), Marquardt and Bruder have studied the elastic and inelastic processes 
in sequential tunneling through two quantum dots in parallel (double-dot interferometer).— 
They have discussed "renormalization effect" in the elastic process and "dephasing" in the 
inelastic process. The former stems from the difference between the ground state of an 
environment in the presence or absence of an electron in the quantum dots. The effect 
reduces the effective tunnel coupling between the quantum dots and leads. This is not the 
case in our resonant tunneling in which an electron virtually stays in a quantum dot, with 
emitting and absorbing phonons, in the elastic process. We do not find the reduction of the 
peak width, whereas the zero-point fluctuation of the phonons decreases the resonant peak 
height. 

The organization of the present paper is as follows. In Sec. II, we present our model for 
the AB interferometer with an embedded quantum dot. The expression of the current is 
derived in terms of the retarded Green function. Section III is devoted to the explanation of 
how to treat the e-ph interaction in the Keldysh Green function formalism. The calculated 
results with optical phonons are shown in Sec. IV, whereas those with acoustic phonons are 
given in Sec. V. The conclusions are presented in Sec. VI. 

II. MODEL AND CALCULATION METHOD 
A. Model 

Our model for the AB interferometer is shown in Fig. ^ There are two paths between 
leads L and R, one path connects the leads through a quantum dot by £l and £r, and the 
other path connects the leads directly by We ltp (reference arm). The phase <p represents the 
magnetic flux in the ring. The bias voltage between the leads is given by eV = hl — jiR, 
where /i£ (//r) is the chemical potential in lead L (R). We fix hl = eV and /i# = 0. Omitting 
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FIG. 1: Model that we adopt for an AB ring with an embedded quantum dot in the presence of 
electron-phonon interaction inside the quantum dot. The phase ip represents the magnetic flux in 
the ring. 

the spin indices, the Hamiltonian for electrons is written as 

H el = H L + H R + H T + H D , (1) 
H L (R) = ^2^kcl {R)k c L{R)k , (2) 

k 

Hd = e d^d, (3) 
H T = 5> L 4 fc d + H.c.) 

k 

k 

+ ^(W^ fc ,c Lfc + H.c), (4) 

kk' 

where cL R ) fe and cl(r)a; denote the creation and annihilation operators of electron with 
momentum k in lead L (R), respectively. We assume a single energy level Eq in the quantum 
dot, with creation and annihilation operators being d\ d. The electron-electron interaction 
is neglected. 

The strength of the tunnel coupling between the quantum dot and leads is characterized 
by the level broadening, T = T L + T R with T L ^ = 2irt\, R -,v, where v is the density of states 
in the leads. We define £ = -k 2 v 2 W 2 for the direct tunnel coupling between the leads. The 
transmission probability through the reference arm is given by T r = 4£/(l + £) 2 - We examine 
the Breit-Wigner resonance without the reference arm by choosing £ = and Fano-resonance 
by 0. 
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We consider the e-ph interaction inside the quantum dot. The Hamiltonian for the 
interaction and that for the phonons are given by 

H e _ ph = ^M q {a q + a ] _ q )Sd, (5) 

H ph = ^uj q a\a q , (6) 

respectively, where a) (a q ) creates (annihilates) a phonon with momentum q. Unless other- 
wise noted, we will be setting h to unity. The coupling coefficient M q is written as 

M q = X q (d\e^ r \d), (7) 

with A g being the amplitude of the e-ph interaction in the two dimensional electron gas. \d) 
represents the envelope function of electrons in the quantum dot. It should be noted that 
the electron-phonon interaction is negligibly small when \q\ > 27r/L, where L is the size of 
the dot, owing to an oscillating factor of (d\e iq ' v \d). Hence we can restrict the summation 
over q to be \q\ < 2n/L in Eqs. (0) and (jOJ. 
The current is expressed as [see Appendix A] 



2e 
7T 

x 



du[f L (u) - f R (uj)] 



T r + v/aT r (l - T r )fcos(^ReG^ d (cj) 

~ [a( 1 - T r cos 2 <p) - T r ] flmG^ dd (o;) J , (8) 

where G r dd (u) is the Fourier transform of the retarded Green function of the quantum dot, 

G r dd (t-t')=-i6(t-t'){{d(t),d\t')}). (9) 

/z,(o;) is the Fermi distribution function in lead L [R]. The level broadening is renor- 

malized asr = r/(l + £). a = ^lYr/V 2 is the asymmetric factor of the quantum dot. 

We consider the tunnel couplings in exactly before taking the e-ph interaction into 
account. In the absence of the interaction, the retarded Green function is given by the 
equation-of- motion method [see Appedendix A], 

Cdd 0) H = -■ ( 10 ) 

uj - e + aT T T cos f + |r 



6 



-4 -2 2 4 

(eo-n L )/r 



FIG. 2: The differential conductance as a function of the dot level £o in the absence of e-ph 
interaction. The phase by the magnetic flux inside the ring is (p = (dotted line), <p = 7r/2 (dashed 
line), and (p = ir (solid line). We fix £ = -k 2 v 2 W 2 = 0.05 (T r = 0.18) and asymmetric factor a = 1. 

The differential conductance dl/dV at zero temperature is expressed as 

^4^, (ID 

1 h e 2 + Y 2 /A V ' 

where e = fi^ — e + jy/aT T T costpM This is an extended Fano form 5 with a complex 
parameter q F , 

q F = yTfiVl ~ T r cos (p + i sin (p) . (12) 

Figure |2] shows the differential conductance as a function of the dot level Eq, with a = 1 
and £ = 0.05 (T r = 0.18). We observe a Fano resonance with dip and peak when (p — 0. 
With increasing the asymmetric resonant-shape changes to symmetric (ip = vr/2), and to 
asymmetric with peak and dip (ip = n). This is in agreement with the experimental results.-^ 



III. TREATMENT OF E-PH INTERACTION 



We consider the e-ph interaction by the perturbation expansion; second-order perturba- 
tion in the case of optical phonons and self-consistent Born approximation in the case of 
acoustic phonons . 17 ! 18 For the comparison, we also examine the Breit-Wigner resonance with 
optical phonons by the canonical transformation method . 2 ° l23 
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A. Keldysh Green function 

To examine the current under a finite bias, we use the Keldysh Green function formalism. 
The Green function of the quantum dot is defined by 

G dd {uj) = I _ (13) 

where the elements are the Fourier transforms of 

G*dd(*-0 = -iT(d(t)d\t')), (14) 

G? dd (t-t') = -if(d(t)S(t')), (15) 

G< d (t-t>) = i(d\t')d(t)), (16) 

GUt-t') = -i{d{t)d\t')), (17) 

respectively. The unperturbed Green function, Gf d (u), is defined in the same way. They 
satisfy the Dyson equation 28 

G dd (co) = G d %) + GgMEHGddM, (18) 

where E(u;) is the self-energy by the e-ph interaction 

S(w = 1 ^ } ) \. (19) 
\e» -E'(u;) ; 

We need the retarded Green function G dd (a;) in Eq. (JBJ) to obtain the transport properties. 
It is expressed as 

G ' dd {uj) = G dd {uj) - G^ d (u). (20) 
G dd (u) follows the Dyson equation of 

G5» = G r d f(u) + G r £\u)^(u)G dd (u), (21) 
where the self-energy is given by 

£ r (u;) = -S<(cj). (22) 

Since G^jj (a;) is given by Eq. (fTUj) . Eq. (|2T|) yields the expression of 

G dd M = 1 - 1 . (23) 

uj — e + jy/aT r T cosy? + ^iT — X r (co>) 

By calculating E*(o>) and S < (a;), we obtain E r (o>) by Eq. (|2*2^1 and hence G^ d (u;) by Eq. 



B. Perturbation expansion 



In the second-order perturbation of e-ph interaction, the elements of the self-energy, Eq. 
(|19|) . are written as 

S » = L E f ^diV - u')D\q^) (24) 

q J 

(7 stands for t, <, >). To obtain E*(o;) and S < (c<j), we need G aa \uj), G dd °' ) (a;) and Green 
functions of phonon, -D'(q, u), D < (q,u). 

G aa °\uj) is calculated using the equation-of-motion method [see Appendix A] as 



G5?V) = G^Mj-fv^sin^^M-^M] 



■ (1 + c)2 [(Fl + FdfLiv) + (I* + ^)/*(o;)]}GdSV), (25) 
where G^u;) is given by Eq. (HHD and G dd °V) = [G dd °V)]*- G* ( d 0) (u;) is given by the 



■■<4alioi..G£ ) ( W ) = Gi 0) (a;) + G< 0) (a;). 

The Fourier transforms of phonon Green functions are 

D\q,u) = -2iri[N q 8(u + uj q ) + N q 8(u - aj q )\ 



+ 77 77, (26) 

LO — LOq + 10 UJ + UJ q — ZO 

D<{q,uo) = -2m[(N q + l)8(uo + uo q ) 

+N g 5(u-u q )]. (27) 

N q denotes the phonon occupation number, l/[exp(u q /kBT) — 1]. We fix the temperature 
at T = in the calculations. 

For the self-consistent Born approximation, G^\u — to') in Eq. (}2*4"j) is replaced by 

G dd(^ - w as 

S'H = ^El M ^| 2 / dw'G^iu-^D^u'), (28) 



and 



S< (^ = ^El M ^| 2 / d^G^a;- 00^(9, wO- (29) 
From the Dyson equation, Eq. (fT8|) . G dd (co>) is expressed as 

G^M = [1 + G dd (u;)E>)]G< (0 V)[l + ^H^dH] + G^E^G^u,). (30) 
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The substitution of Eqs. flU and (j2SJ) into Eq. flU)) yields 



[w - £ + Iv^rcos^ - ReE r (u;)] 2 + [|r - ImS r (w)] : 



x ^-r^aT lS m<p[f L (u)-f R (u)\ 

G^ d (u;) is obtained by Eq. flUJ), using G< d {u) in Eq. (JSTJ) and G r dd {uj) in Eq. (J23D- We 
determine £ r (u;) and S < (w) by solving the equations self-consistently. 



C. Canonical transformation 

For the comparison with the perturbative method, we examine the canonical transforma- 
tion method . 20 ! 23 In this method, we decouple the e-ph interaction in Eq. (jSj) by a canonical 
transformation o£2& 

H = e s He- s , (32) 
s = d^d^^ — -{a)_ q — a q ). (33) 

Q 9 

By the transformation, the Hamiltonians for electrons in the quantum dot, e-ph interaction, 
and phonons are rewritten as 

u q a\a q . (34) 

Q Q 

The tunneling Hamiltonian, Eq. (J3J), is transformed to 

S T = ]>> L c[ k J+H.c.) 

k 

+ I](^RCRk^+H.c.) 

k 

+ ^(H/e^4 k ,c Lk + H.c), (35) 

k,k' 

where d and dft are renormalized operators of electrons, 

d = e s de~ s = dX, (36) 
10 



d ] = e s $e 



s 



(37) 



with 



X = exp -E^K 




(38) 



After the decoupling of electron and phonon operators in Eq. (J34|) . we take into account the 
e-ph interaction exactly, whereas we adopt an approximation for H^: Electron operators, d 
and d\ are replaced by d and d^, disregarding the e-ph interaction in the tunnel couplings. 
The retraded Green function G dd (t — t') is given in Appendix C. 

IV. OPTICAL PHONONS 

We begin with the case of longitudinal optical phonons. To illustrate the decoherence 
— zero-point fluctuation effect in the elastic process and dephasing effect in the inelastic 
process — , we mainly discuss the effects on the Breit-Wigner resonance in this section. 

The amplitude of the e-ph interaction is described by the Frohlich-coupling 28 



where m* is the effective mass of electrons. e(oo) and e(0) are the dielectric constants at 
high and low frequencies, respectively. V is the normalization volume. As a good approxi- 
mation, we assume that all the phonons have the same energy u>o (Einstein model) since the 
wavenumbers are limited to a small range of \q\ < 2-7r/L, when L is in the submicron scale. 
We define the coupling strength of the e-ph interaction as 



A. Breit-Wigner resonance 

We examine the Breit-Wigner resonance through a quantum dot without the reference 
arm (£ = 0, T r = 0). We treat the e-ph interaction by the second-order perturbation. In Fig. 
El we plot the differential conductance, dl/dV, as a function of the energy level e in the 




(39) 



(40) 




(41) 
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quantum dot. The bias-voltage is eV = 0.5T (eV < ujq = 2T; dashed line) and 3r [eV > u>o; 
solid line). A dotted line indicates the conductance in the absence of the e-ph interaction, 
where the conductance is independent of the bias. A subpeak at e — /il ~ — u> is seen only 
when eV > Uq, which corresponds to the emission of a phonon. 

When eV < uq, a real process of phonon emission cannot take place at T = 0. Even in 
this case, the peak height is suppressed by the e-ph interaction. This is due to the elastic 
process in which electrons emit and absorb phonons virtually. It disturbs the coherent 
transport through the quantum dot and hence diminishes the resonant peak height. When 
eV > ujq, both elastic and inelastic processes are possible. The latter process makes the 
subpeak of the phonon emission, whereas the former process mainly suppresses the main 
peak, as discussed in the next section. 

The position of the main peak is shifted by the e-ph interaction, regardless of eV < u>q 
or eV > ujq. This is ascribable to the above-mentioned elastic process. The zero-point 
fluctuation of phonons results in the renormalization of the energy level in the quantum dot, 
in addition to the suppression of the peak height. 

The width of the main peak is influenced by the inelastic process of phonon emission: 
it does not change when eV < uj , whereas it is broadened when eV > uj although it is 
not clearly seen in Fig. EJ To discuss the width, we examine the retarded Green function 
in Eq. ()23|) . The peak width is determined by the imaginary part of the self-energy S r (w) 
in its denominator. When eV < ujq, Im£ r (co>) = in the range of integration, < u < eV, 
in Eq. (jHJ). Hence the peak width is not changed. We conclude that the peak width is not 
influenced by the elastic process of virtual emission and absorption of phonons in the resonant 
tunneling. This is in contrast to the sequential tunneling where the renormalization of the 
phonon state reduces the effective tunnel coupling between the quantum dot and leads.«2A 
When eV > uo, on the other hand, ImE r (ct;) ^ 0. At uj = eV, 

c 2 r 

ImZ r (eV) = - — - ^ — , (42) 

which mainly determines the linewidth of the differential conductance. Since ImE r (eV) < 0, 
the linewidth is broadened by the e-ph interaction. The real emission of phonons decreases 
the life-time of the dot state and, as a result, increases the peak width. 
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FIG. 3: The differential conductance as a function of the dot level Eo> hi the case of Breit-Wigner 
resonance (£ = 0, T r = 0) with optical phonons. The asymmetric factor is a. = 1. The phonon 
energy is wo = 2r, whereas the strength of e-ph coupling is £ = 0.8r. The bias voltage is eV = 0.5T 
(dashed line) and 3r (solid line). Dotted line indicates the conductance in the absence of e-ph 
interaction. 

B. Elastic and inelastic processes 

To investigate the elastic and the inelastic processes more precisely, we perform the 
calculation considering the elastic process only. Among the four elements in the self-energy 
E(kj), Eq. (jl9|) . the diagonal elements, E* and £*, describe the elastic process while the 
off-diagonal elements, and E > , represent the inelastic process, as discussed in Appendix 
B. Excluding the latter, the retarded Green function is written as 



In Fig. |U we present the calculated results using Eq. (J4*3~j) by dashed lines. Solid lines 
indicate the results in the presence of both elastic and inelastic processes. When eV < ujq 
[Fig. 13(a)], both the results are identical to each other since only the elastic process exists. 
When eV > ujq [Fig. 13(b)], the inelastic process makes some differences. The height of the 
dashed line and solid line are almost the same, which indicates that the suppression of the 
main peak is mostly caused by the elastic process. Note that the elastic process also results 
in a subpeak at Eq — /xl ~ when eV > ujq. This is due to the formation of so-called 





(43) 
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FIG. 4: Differential conductance as a function of the dot level Eq in the case of Breit-Wigner 
resonance (£ = 0, T r = 0) with optical phonons. The asymmetric factor is a = 1. The phonon 
energy is u>o = 2r, whereas the strength of e-ph coupling is £ = 0.81~\ The bias voltage is (a) 
eV = 0.5T and (b) 3I\ The dashed line includes only the elastic process of e-ph interaction. The 
solid line includes both the elastic and inelastic processes. Note that dashed and solid lines are 
overlapped by each other in (a). 

polaron in which an electron is coherently coupled with phonons. The resonant tunneling 
through a polaron level makes the subpeak. The polaron has a finite lifetime by the real 
emission of phonon. Hence the inelastic process broadens the subpeak as well as the main 
peak. 

C. Fano resonance 

Figure shows the Fano resonance with optical phonons, in the presence of reference 
arm (x = 0.05, T r = 0.18). The magnetic phase is (a) ip = 0, (b) ip — ir/2, and (c) ip = n. 
The differential conductance, dl/dV, is plotted as a function of the energy level e in the 
quantum dot, when eV < u (dashed line) or eV > u> (solid line). The dotted line indicates 
the non-interacting case. 

The influence of the e-ph interaction on the Fano resonance is qualitatively the same 
as that on the Breit-Wigner resonance. When eV < ujq, only the elastic process takes 
place. Then we do not observe a subpeak of phonon emission. The zero-point fluctuation 
of phonons reduces the resonant amplitude; both the resonant peak and dip are diminished. 
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FIG. 5: Differential conductance as a function of the dot level Eq in the Fano resonance (£ = 0.05, 
T r = 0.18). The asymmetric factor is a = 1. The phonon energy is loq = 2r, whereas the strength 
of e-ph coupling is ( = 0.8T. The magnetic phase is (a) (p = 0, (b) (p = n/2, and (c) ip = tt. The 
bias voltage is eV = 0.5T (dashed lines) and 3r (solid lines). The dotted lines indicate the case in 
the absence of e-ph interaction. 
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The resonant width is not changed by the elastic process. When eV > ojq, a subpeak appears 
by the real emission of phonon. The inelastic process broadens the width of main resonance. 

The shape of the Fano resonance changes with magnetic phase (p in the same way as in 
Fig. |21 We find that the renormalization of the phase is negligibly small. 



D. Calculation by canonical transformation method 

In this subsection, we consider the e-ph interaction by the canonical transformation 
method. By the method, the retarded Green function is easily obtained a o 23 ' 28 



oo 



G r dd H =e^ N ^ J2 JpkjN q (N q + l)\ 



l=— oo 

e lu /3/2 



uj — (e — A) — Iujq + \\faT r Y cos<^ + ~r 



(44) 



where Ji(x) are the Bessel functions, k = ((/uo) 2 and energy shift is A = C 2 /W At zero 
temperature, G dd (cj) is expressed as 

CT^tj) = e~ k V k ' l \ ^_ ^. (45) 

w - (e - A) - lu + j VaT r T cos (p + |T 

The substitution of Eq. (J4~5j) into Eq. (JEJ) yeilds the differential conductance. In the case 
of Breit-Wigner resonance (£ = 0, T r = 0), 

, TlJTr 2e 2 A r 2 e- k k l /n 

"I* - X £ T ^-fa-A)-^ + P/4 - (46) 



In the case of Fano resonance, 



!=— oo i ' 

where e = yU L — (e — A) — Zc<j + -v^aTIr cos<^. 

The calculated results are shown in Fig. El (a) Breit-Wigner resonance and (b) Fano 
resonance. In the case of Breit-Wigner resonance, the main peak (I = 0) is given by the 
Lorentzian form with center at eq — //l = A, indicating the renormalization of the energy 
level. The peak height is reduced by the factor of e~ k , whereas the width is fixed at T/2. A 
subpeak is seen at e Q — /il = A — u , corresponding to the emission of a phonon. Hence the 
polaron formation can be described by this calculation method, similarly to the Tien-Gordon 
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FIG. 6: Differential conductance as a function of the dot level Eq with optical phonons, calculated 
by the canonical transformation method. The asymmetric factor is a = 1. The phonon energy is 
cjo = 2r, whereas the strength of e-ph coupling is £ = 0.8r. Dotted line indicates the differential 
conductance in the absence of e-ph interaction, (a) Case of resonant tunneling through a quantum 
dot (£ = 0, T r = 0). (b) The case of Fano resonance (£ = 0.05, T r = 0.18). 

theory.— The calculated results of Fano resonance can be understood in the same way as 
those of the Breit-Wigner resonance. 

We find some shortcomings of this method, (i) The bias-voltage dependence of the 
transport cannot be obtained. The dl/dV curve is independent of the bias in Eqs. ()46|) 
and ([47)1 . Hence the inelastic process of phonon emission takes place even when eV < ujq. 
The renormalization of the energy level by the e-ph interaction does not depend on the bias 
voltage, (ii) The resonant width of the main peak is not changed by the phonon emission. 
Hence the finite lifetime by the inelastic process is not considered properly. These are due to 
the approximation in the tunnel Hamiltonian. By the replacement of the electron operator 
d by d in [Eq. (135(1 ]. we neglect the change of the phonon states accompanied by the 
tunnel processes. 

V. ACOUSTIC PHONONS 

Now we discuss the e-ph interaction with acoustic phonons. The amplitude of the e-ph 
interaction in Eq. J7J) is described by the piezoelectric coupling , 22 i 31 




V\q 



(48) 
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FIG. 7: Differential conductance as a function of the dot level £q in the Breit-Wigner resonance 
with acoustic phonons. The asymmetric factor is a = 1. The coupling constant is g = 0.05, 
whereas the dot size is L = c s /(2T). The bias voltage is eV = 0.5T (dashed line) and 217 (solid 
line). The dotted line indicates the conductance in the absence of e-ph interaction. 

where g is a coupling constant and c s is the sound velocity. The dispersion relation of the 
phonons is given by 



in Eq. (J7J). For the calculation of e-ph interaction, we adopt the self-consistent Born approx- 
imation. 

A. Breit-Wigner resonance 

We begin with the Breit-Wigner resonance in the absence of reference arm (£ = 0, T r = 0). 
We set g = 0.05 and c s = 5000m/s, which are the values in GaAs . 19131 We fix L = c s / (2T): e. 
g. when T = O.lmeV, the dot size is L = 0.015/im (smaller than in Ref. [3]). The differential 
conductance is plotted in Fig. as a function of the energy level Eq in the quantum dot. 
The bias voltage is eV = 0.5T (dashed line) and 2T (solid line). In contrast to the case of 
optical phonons (Sec. IV), the inelastic process always exists, which is accompanied by the 



LO q = C s \q\. 



(49) 



For the convenience of the calculation, we assume that | (d\e iq ' r \d) | 2 is given by 




(50) 
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FIG. 8: Differential conductance as a function of the dot level £q in the Fano resonance with 
acoustic phonons. We set £ = 0.05 (T r = 0.18). The asymmetric factor is a = 1. The coupling 
constant is <? = 0.05. (a)The dot size is fixed at L = c s /(2T). The bias voltage is eV = 0.5T 
(dashed line) and 2r (solid line). The dotted line indicates the conductance in the absence of 
e-ph interaction. (b)We change the dot size as L = c s /T (dotted line), c s /(2T) (dashed line), and 
c s /(5r) (solid line). The bias voltage is eV = 2T. 

real emission of phonons. A subpeak structure is not observed since the acoustic phonon 
has a continuous spectrum. 

We find the following effects of the e-ph interaction, (i) The resonant peak height de- 
creases with increasing the bias voltage. As discussed in the previous section, this should 
be mainly ascribable to the elastic process, (ii) The peak width is more broadened with 
the bias voltage, which is due to the inelastic process. The emission of phonons results in 
the finite lifetime and hence an increase in the resonant width, (iii) The form of the Breit- 
Wigner resonance becomes asymmetric. The resonant peak is more suppressed on the side 
of Eq < than on the side of Eq > jii,. This may be explained as follows. When e < fi^, 
phonons can be emitted in the electron tunneling from lead L to the quantum dot, or in the 
tunneling from the dot to lead R. When e q > fi^, phonons can be emitted only in the latter. 
Hence the decoherence works differently on both sides of the resonance. 

B. Fano resonance 



Figure IHfa) shows the Fano resonance with acoustic phonons, in the presence of reference 
arm (£ = 0.05, T r = 0.18). The bias voltage is eV = 0.5T (dashed line) and 2T (solid line). 
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A subpeak structure is not seen as in the Breit-Wigner resonance with acoustic phonons. 

The e-ph interaction decreases the resonant amplitude, whereas it increases the resonant 
width, reflecting the real emission of phonons. These effects are more prominent with larger 
bias. Particularly, the resonant dip becomes almost invisible with high bias and as a result, 
the asymmetric shape of the Fano resonance grows like a symmetric one. Although this 
is in qualitative accordance with the experimental results,— the resonant peak is still in an 
asymmetric form even in a smaller size of quantum dots (with stronger e-ph interaction, 
as discussed below) than in the experiment. To explain the experimental results quantita- 
tively, we would have to consider the e-ph interaction in the reference arm as well as inside 
the quantum dot. The other dephasing effects, e.g. spin-flip of a localized electron in the 
quantum dot, might be also important. 

We change the size of the quantum dot, L, in Fig. |H[b): L = c s /T (dotted line), c s /(2T) 
(dashed line), and c s /(5T) (solid line). The bias voltage is fixed at eV = 2T. With an 
decrease in the dot size, the matrix element (d\e iq ' r \d) in the coupling coefficient M q increases 
[see Eq. (jHUj) ] . In consequence, the effect of e-ph interaction is enhanced. The figure indicates 
smaller amplitudes and wider widths of the resonance in smaller quantum dots. 

VI. CONCULSIONS 

We have investigated the Breit-Wigner resonance and Fano resonance in quantum dots 
with e-ph interaction. Using Keldysh Green function method and perturbation with respect 
to the e-ph interaction, we have obtained the bias- voltage (V) dependence of the decoherence. 

With optical phonons of energy uq, we distinguish elastic and inelastic processes. When 
eV < Uq, only the elastic process takes place in which electrons emit and absorb phonons 
virtually. The process suppresses the resonant amplitude, whereas it does not change the 
resonant width. We do not observe the renormalization of the tunnel coupling which has 
been discussed in the sequential tunneling.— When eV > u , the inelastic process is possible, 
which broadens the resonant width. The process also makes a subpeak, as a function of the 
energy level in the quantum dot, corresponding to the real emission of phonons. We cannot 
obtain any bias- voltage dependence of the decoherence when we consider the e-ph interaction 
by the canonical transformation method, neglecting its effect on the tunnel coupling. 

With acoustic phonons, the suppression of resonant amplitude and broadening of the 
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resonant width are more prominent with an increase in the bias. In the Fano resonance, 
the resonant dip becomes almost invisible. The asymmetric resonant shape becomes like 
a symmetric one, in qualitative accordance with the experimental results.— However, the 
model with e-ph interaction in the quantum dot only is not sufficient for the quantitative 
explanation of the experiment. 

In the present paper, we use the term of dephasing for the real emission of phonons in 
the inelastic process. This is because the elastic process does not result in the dephasing, as 
discussed in Ref. |3 in the presence of reference arm, one might think that the dephasing 
effect can be directly evaluated by calculating the amplitude of AB oscillation. But it is not 
correct. The e-ph interaction leads to the elastic and inelastic processes in the quantum dot 
on one hand, it changes the ratio of the current through the reference arm to that through 
the quantum dot on the other hand. Hence the amplitude of the AB oscillation cannot 
measure the dephasing effect correctly. 

The authors gratefully acknowledge discussions with F. Marquardt. This work was par- 
tially supported by a Grant-in- Aid for Scientific Research in Priority Areas "Semiconductor 
Nanospintronics" (No. 14076216) of the Ministry of Education, Culture, Sports, Science and 
Technology, Japan. 

APPENDIX A: EXPRESSION OF CURRENT 

The current from lead L to the quantum dot is given by the time derivative of the 
occupation number in lead L, 



e(A L ) 



ie{[H,N L ]), 



(Al) 



with Al = X]fc c Lfe c Lfc- A straightforward calculation yields 



h = ia 





(A2) 



Here, we have introduced lesser Green functions, G^ L 
G^ k , Lk (t,t') = i(c Lfe (£')cRfc/(t)). In the energy representation, 



(MO 



i(cl k (t')d(t)) and 



J L (w) = -7- Re / du [t h + We~* £ G< fe , hk (u)] 



(A3) 
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The current from lead R to the dot, Ir, is obtained in the same way. In the case of non- 
interaction leads, we can express the current using the Green function of the quantum dot 
only, as shown in the following. 

Let us consider G^ Lk (u)) in Eq. (jA3|) . Using equation-of- motion method and analytic 
continuation rules,— G dLk (t — t') is written as 

G< Lk (t -t>) = J d tl [t L G r dd (t - fi)0&(ti - 
+ ttG^t-^g^ih-t') 
+ We^G^t-h^h-?) 

+ We^G< Rk (t- tl )gt k ( tl -t% (A4) 

where G d , Rk (t - t') = -i9(t - t')({d(t), G< Rk (t - t') = i( CKk (t')d(t)), fa(t - if) = 

i9(—t + t') ({cLk(t) , c^ Lk (t')}) , and g^ k (t — t') = «(cl A (^') c lA;(0) ■ The Fourier transformation 
yeilds 

53 G d,LkM = iKvt h G^ A {uj) + 2mvf h (uj)t h G 7 Ad {uj) 

k 

+ mvWe^Gi,^) +MMW^(7 d », (A5) 

k k 

where J^.^^) = inu and J2 k g^ k (uj) = 2mvf L (uj). A similar expression of Y, k G dRk (u) 
is obtained in the same way. The expression together with Eq. (jA5|) yeilds 



k ^ 



+ Y--[t U K)fHK)(u) +i™We^n R{h) f R{h) {u)]G dA {u) 
2t7ruWe^f L(R) (uj) ^^ r 

+ [Tt 2^ G d,R(L)fcl^J 

7—7 2^ G d,L(R)fe (w) • ( A6) 

^dL(R)fc( w )' which appears on the right-hand side of Eq. (|A6|) . is derived from the 
equation-of-motion method as 

^d,L(R)fe(^ ~ = J dti[tt( R )G dd (t - tl)fl , £(R)fe(*l ~~ 
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where <7£(R)jfe(^ — = — ^6>(t ~~ ({ c L(R)fc(£)j c L(R)fc(^')})- From the Fourier transform, we 
obtain 

E G d,L ( R )fc = t^tH^lcr) - ^^e ±! \ (L) ]^ d ( W ). (A7) 

Then we express G dL , R s k (ui) in terms of G r dd (u) and G^ d (u;). 

G Rk , Lk (uj) in Eq. ()A3|) can be expressed using G dd (u) and G dd (u), by the same technique. 
As a result, the current /l(^) [and also /r(u;)] is rewritten in terms of G dd (u) and G dd (u) 
only. Finally we delete G dd (u) by taking a linear combination of them, I = x/l — (1 ~ ^Kr, 
and obtain Eq. (JSJ). 

The expression of unperturbed Green functions, Eqs. (fTUj) and (J23j) . are obtained also by 
the equation-of-motion method: 

G dd 0) (t-t') = j d tl [g r d ( tl - t>) 

+tLG r d %(t- tl )f d ( tl -t') 

+tRG r $ k (t-t 1 )g r d (t 1 -t% (A8) 

GT(t-t') = J dhig^h-t') 

+tLG r d %(t- tl )g<( tl -t') 

+^G d <i°i(t-tikd(ti-t') 
+tnG r d %(t-t 1 )g<(t 1 -t') 

+tRG^ k (t-t 1 )g a d (t 1 -t% (A9) 

where g d (ti — £') and <? d (ti — £') are the Green functions of the non-interacting electrons 
in an isolated dot, 

<7d ( °V) = „ * ■ (MO) 

3d( w ) = -2nirid5(uj - £ )- (All) 
rid is the distribution function of electrons in the dot. By the Fourier transforms, we obtain 
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FIG. 9: Complex-time contour on which nonequilibrium- Green function theory is constructed. 

fit t I oo fi^ t , ==T) oo 

X' T 

G ' ' | > ) 00 G I , ' > ) 00 

FIG. 10: Diagrams of the coutour-ordered Green function G(t,t'). Four diagrams correspond to 
the real-time Green functions, Gr, G l , G > , and G < , respectively. 

+i R <S ) fe H< ?d <(u,) + tKG<^ k (u)g a d (u). (A13) 

By expressing gJJ^w), GjS (<*>), ^b*M> and G d,klM in terms of G ddV) and G dd (0) H, 
we obtain Eqs. flE) and Eq. (f25|) . 



APPENDIX B: PERTURBATION EXPANSION 



We discuss elastic and inelastic processes by the e-ph interaction, using the coutour- 
ordered Green function. The contour-ordered Green function of the quantum dot is defined 

as 

G dd (r,r') = -i{T c {Sd{r)S{r')}), (Bl) 

where 



S = Tr 



c\ exp 



c 



(B2) 



is the contour-ordered S matrix and d(r) = e lHr de~ lHr . As shown in Fig. the contour 
path C goes from t = — oo to t = oo and goes back to t — —ooM The time labels r and r' 
are located on the countour path. The contour-ordering operator Tq orders the operators 
following it in contour sence: Operators with time labels later on the contour are moved left 
of operators of earlier time labels. 
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T 2 \ /T1 

L t '-■ ■ -, > ) 00 L < t - ■ -" > ) 00 

FIG. 11: Diagrams of the coutour-ordered Green function with self-energies. Solid lines indicate 
the contour-ordered Green function of the dot electron, whereas dashed lines indicate the contour- 
ordered Green function of phonons. Four diagrams correspond to the real-time Green functions 
with self energies of £*, £*, S < , and T, > , respectively. 

The contour-ordered Green function, G(t,t'), is related to the real-time nonequilibrium 
Green functions, as shown in Fig. G l (t — t') when r and r' are on the upper branch, 
G l (t — t') when r and t' are on the lower branch, and G > (t — t') [G < (t — t')] when t' is on 
the upper [lower] branch and r is on the lower [upper] branch. 

The contour-ordered Green function obeys Wick's theorems in the contour-time order. 
In the second-order perturbation with respect to the e-ph interaction, 

G dd (r,r') = G° dd (r,r')+ [ dr x [ dr 2 G° dd (r, n)E(ri, r 2 )G dd (r 2) r'), (B3) 

Jc Jc 

= M q G° dd (n, T 2 )D(q, r 1; r 2 ), (B4) 
q 

where 

GUn, t 2 ) = -t(T c {d(r)d\r')}} eh (B5) 
D(q, n, r 2 ) = -i^cQ^Q^))^ (B6) 

with 

Q(r)=a q (r) + al q (r). (B7) 

The bracket ( ) e i denotes the average over the states of dot electron in the absence of 
e-ph interaction and ( ) p h denotes the average over those of phonons in the equilibrium 
state. Figure ^2 shows four types of diagrams for the self-energies. Solid lines indicate 
the contour-ordered Green function of the dot electron, whereas dashed lines indicate the 
contour-ordered Green function of phonons. The diagrams can be divided into two types, 
one corresponds to the elastic process and the other corresponds to the inelastic process. 
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When T\ and r 2 are on the same branch of contour path C, electrons emit (absorb) and 
absorb (emit) phonons. Therefore, the phonon state at t = 00 is the same as at t — —00. 
This is the elastic process in which the energy of electrons conserves. On the other hand, 
when 7~i and r 2 are on the different branches of the path, electrons emit or absorb phonons 
before t = 00. This process indicates the real emission or absorption of phonons and hence 
corresponds to the inelastic process. Thus, the processes are inelastic. 

In pratical calculations, we have adopted the real-time nonequilibrium Green functions. 
These four diagrams correspond to self-energies of real-time nonequilibrium Green fuctions 
£*, E*, S < , and S > , as labeled in Fig. ^2 Note that our definition of elastic and inelastic 
processes differs from that in Refs. 19 and 20. 

APPENDIX C: CANONICAL TRANSFORMATION 

By the canonical transformation of 20,23 

H = e s He~ s (CI) 
and approximation of the replacement of by H^, we can divide the Hamiltonian into an 

Ha = [eo ~ J2( M i/ u ^ d + H T + H L + H R (C2) 

Q 

and phonon part = H^. The retarded Green function is decoupled as 

Gdd(*~0 = -i9(t-t')({d(t)J^t')})4X(t)X\t')) ph , (C3) 

where 

d(t) = e i£ ^de- in ^ (C4) 

and 

X(t) = e iS v ht Xe- iS ? ht . (C5) 

The bracket ( ) e i denotes the average over the states of dot electron in the absence of e-ph 
interaction and ( ) p h denotes the average over those of phonons in the equilibrium state. 

The electron part is calculated by the equation-of-motion method [similar to Eq. (|A8j) ] 
and the phonon part is calculated using Feyman disentangling.^ We obtain the retarded 
Green function of the quantum dot 

G r dd (t-t>) = G^ l \t-t')F(t-t'), (C6) 
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where 



Gdd Gl) (* - = ~W(t ~ exp{-t[e - ^(M 2 /o; q ) - - JdT x Y cos <p]t + f/2t} 



and 




(C7) 



(C8) 



(C9) 



The Fourier transform of G dd (t — t') yeilds Eq. 
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